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We report the current status of the observations of supernova remnants (SNRs) with the Large 
Area Telescope (LAT) aboard the Fermi Gamma-ray Space Telescope, focusing on middle-aged 
SNRs that appear to be interacting with molecular clouds. Observations with the Fermi LAT in an 
energy range from 0.2 GeV to ~ 100 GeV have unveiled the presence of luminous GeV gamma- 
ray emission in middle-aged SNRs, providing a new insight into the shock-acceleration theory 
and the origin of galactic cosmic rays. The middle-aged SNRs detected by the Fermi LAT are 
generally much brighter in GeV than in TeV in terms of energy flux, which emphasizes the im- 
portance of the GeV observations. Spectral steepening in the Fermi-LAT band is commonly found 
for the GeV-luminous SNRs. Remarkably, most (if not all) of the GeV-luminous SNRs are known 
to be interacting with molecular clouds, and they are also the strong sources of synchrotron radio 
emission. We discuss possible scenarios to explain the enhanced GeV gamma-ray emission in the 
cloud-interacting SNRs. Particular emphasis is placed on a scenario in which shock-accelerated 
cosmic-rays are adiabatically compressed and energized as result of radiative cooling behind the 
cloud shock. 
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1. Introduction 

Diffuse shock acceleration §\§ operating at a collisionless shock produces nonthermal parti- 
cles in a supernova remnant (SNR). Remarkable progress in the observations of shock-accelerated 
particles was made in the last decade, especially in the X-ray and TeV y-ray bands, which has 
strengthened the conjecture that shock acceleration in SNRs is responsible for the bulk of the galac- 
tic cosmic rays (GCRs). Synchrotron X-ray emission produced by multi-TeV electrons has been 
detected from most of young historical SNRs 0. The X-ray measurements by the Chandra X-ray 
Observatory with its arcsecond resolution imaging have shown that particle acceleration to TeV 
energies at strong shocks in young SNRs accompanies significant amplification of interstellar mag- 
netic fields [g]. Direct evidence for the production of multi-TeV particles in shell-type SNRs came 
from the TeV y-ray observations made with ground-based Cherenkov telescopes such as H.E.S.S. 
[Q], though the emission mechanism ("hadronic vs leptonic") remains unsettled. Gamma-ray ob- 
servations at the GeV band have been expected to provide complementary information on particle 
acceleration at the expanding shock waves of galactic SNRs. 

In the GeV y-ray band, possible associations of y-ray sources with five radio-bright shell-type 
SNRs were reported by Ref. [g| based on the data taken in 1991-1994 by the EGRET instru- 
ment onboard the Compton Gamma Ray Observatory (see also Ref. The two strongest y-ray 
sources that are possibly associated with SNRs are the EGRET unidentified sources in the direc- 
tions of y Cygni and IC 443, both of which are known to be interacting with adjacent molecular 
clouds. If the y-ray excesses are due to cosmic-ray interactions with a dense cloud, the dominant 
emission processes are either bremsstrahlung by relativistic electrons or 7T°-decay y-ray s resulting 
from proton-proton hadronic collisions. However, the EGRET point-spread function is much wider 
than the spatial extent of most of the SNRs with possible y-ray associations, which made it difficult 
to preclude the possibility that the y-ray sources are due to pulsars 1 . 

The Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope [Js[] , which 
was launched in June 2008, has started to explore the GeV y-ray sky with much improved resolving 
power and unprecedented sensitivity. The observations with the Fermi-LAI provide new methods 
of addressing unsolved problems in shock acceleration theory and the origin of the GCRs. Us- 
ing the LAT data from the first ~ 1 yr of observations, the Fermi LAT Collaboration has reported 



the detections of y-ray emission from several SNRs [g]- 14]. Middle-aged SNRs interacting with 
molecular clouds (MCs) constitute the dominant class of y-ray-luminous SNRs. Here we summa- 
rize the results for such MC -interacting SNRs obtained with the LAT and discuss the mechanisms 
by which y-ray emission can be enhanced in such systems. 

2. Fermi-LAT Detections of SNRs Interacting with Molecular Clouds 

The Fermi LAT Collaboration has so far reported the discoveries of spatially extended y-ray 
sources coincident with four remnants: W51C @, W44 [|TTJ, IC 443 [jTJ], and W28 [JT3|] . They 



are middle-aged (~ 10 4 yr) remnants that appear to be interacting with molecular clouds. Figure 1 
shows a collection of the 2-10 GeV y-ray count maps obtained with 2.5 yr data of the LAT. For 



1 One of the two most significant EGRET sources (a source in the y Cygni SNR) is indeed found to be a y-ray pulsar 
with the Fermi LAT [|] . 
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Figure 1: Fermi LAT count maps in 2-10 GeV of the MC-interacting SNRs with extended y-ray emission: 
(a) W51C; (b) W44; (c) IC443; and (d) W28. The LAT count maps (in the celestial coordinates) constructed 
from front-converted events in the four panels are smoothed by a Gaussian kernel of 0. 15 deg. The intensity 
is in units of counts per pixel with a pixel size of 0.05 x 0.05 deg 2 , and shown in a common range of 1 .0-5.5 
counts per pixel. Superposed are the contours from the VLA radio maps (see the text). A black ellipse in 
panel (a) represents the location of shocked CO clumps JIB]]. The black crosses in panel (c) are the locations 
of shocked molecular clumps from which OH maser is detected [[To]]. The position of a pulsar wind nebula 
is also marked in panel (c). 



better spatial resolution, only front-converted events are used. Even without subtracting the intense 
Galactic diffuse emission, significant y-ray emission clearly stands out in excess of the background 
diffuse emission. In each panel, superimposed is a radio continuum map synthesized with the 
Very Large Array (VLA) at a frequency of: (a) 1.4 GHz [0; (b) 1.4 GHz (The FIRST Survey); 
(c) 330 MHz and (d) 330 MHz. The spatial extent of the y-ray emission estimated through 
the maximum likelihood technique was generally found to be consistent with the size of the radio 
remnant. Given the size of a point-spread function of the LAT in this energy band, dedicated 
investigations are necessary to establish possible sub-structures within these remnants, which will 
be reported elsewhere. 

The radio and y-ray properties of the four SNRs are similar in many aspects. The synchrotron 
radio emission has a large flux of 160-310 Jy at 1 GHz with flat spectral index of a ~ 0.3-0.4 
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Figure 2: Gamma-ray spectra of shell-type SNRs measured with the Fermi LAT: W51C (red filled circles 
[[§); W44 (red open circles [|TJ); IC 443 (red filled rectangles [|l|]); W28 (red open rectangles 
Cassiopeia A (black filled circles [fffl). 



P18|]. The GeV y-ray spectrum commonly exhibits a spectral break at around 1-10 GeV, and the 
luminosity ranges L y = 10 35-36 erg s . In Fig. |2[ the /-ray spectra of the MC-interacting SNRs 
measured with the LAT are shown in the so-called vL v form along with the spectrum of Cas A. 
The distances are taken from the nominal values described in Refs. [|9]-[l3|]. Note that the y- 
ray luminosities of W51C and W44 are more than an order of magnitude larger than the y-ray 
luminosity of Cas A. 

The LAT-detected SNRs are generally radio-bright objects, which suggests some physical link 
between the synchrotron radiation and the GeV y-rays. Indeed, the median value of the radio 
surface brightness 2 at 1 GHz of the Galactic SNRs is £ ~ 2 x 10~ 9 erg cm~ 2 s _1 sr _1 [18], which 



is an order of magnitude smaller than that of the LAT-detected SNRs. Moreover, there found a 
possible correlation between the radio and y-ray fluxes from MC-SNRs. Figure [3] (left) shows a 
plot of the radio flux vs y-ray flux for the current list of MC-SNRs that are detected by the LAT 
[§-|l4], |i9[ |. The radio and y-ray energy fluxes represent the total fluxes integrated over a remnant. 
In Fig. H (right), a comparison between the radio and y-ray intensities is made in terms of mean 
surface brightness, where the total flux is divided by the solid angle subtended by each remnant. 
The mean surface brightness, which is distance-independent, has a large dynamic range over more 
than two orders of magnitude. 



2 Usually, the radio surface brightness is represented as L v = f v /SI in units of W m 2 Hz 1 sr 1 . Here we adopt 
Z = vS v = vf v /CI in units of erg cm~ 2 s~' sr~' . 
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Figure 3: (Left) Radio flux (synchrotron) vs GeV y-ray flux for MC-interacting SNRs. The y-ray energy flux 
integrated over 0. 1-100 GeV and the radio flux, vf v at 1 GHz, are shown. (Right) Mean surface brightness 
of the synchrotron radio emission and GeV y-ray emission. The flux-flux plot is converted into this form 
using the solid angles of the radio remnants. 



3. Discussion 

Since the interaction with a molecular cloud clearly plays a key role in enhancing the y-ray 
emission, bremsstrahlung by relativistic electrons or 7T°-decay y-rays produced by high-energy 
protons are the two plausible dominant channels of the y-radiation. The observed large luminosity 
of the GeV y-ray emission precludes the inverse-Compton scattering off the CMB and interstellar 
radiation fields as the main emission mechanism. For instance, the luminosity of SNR W51C 
requires an unacceptably large amount of the total CR electron energy of W e ~ 1 x 10 51 ergs [g]. 
Generally the high-energy protons are thought to be the predominant component of accelerated 
particles in SNRs so that the 7T°-decay y-ray emission is stronger than the electron bremsstrahlung 



(see however Ref. [20]) 



There are two different types of scenarios to explain the GeV y-ray emission arising from 
the MC-SNR systems. The Runaway CR model considers y-ray emission from molecular clouds 
illuminated by runaway CRs that have escaped from their accelerators [|T]-^3|. CR particles 
released by a remnant diffuse into interstellar space with an energy-dependent diffusion coeffi- 
cient, and collide with massive clouds producing y-ray photons via 7T°-decay (by CR hadrons) and 
bremsstrahlung (by CR electrons) channels. The y-ray wide-band spectra from the GeV to TeV 
energies can markedly vary from one object to another depending on the relative locations of the 
accelerator and the target cloud as well as the time history of the accelerator. For example, a hard 
GeV y-ray source [|^] which is also bright in TeV [24] outside the southern boundary of SNR 
W28 may be explained by runaway CRs. Future analysis of spatial distributions of GeV y-rays 
combined with CO observations would make it possible to discern the y-ray sources that are at- 
tributable to runaway CRs. Then one can learn about how CRs are released into interstellar space 
and how they propagate in the vicinity of the SNR where self-generated Alfven waves may change 



the local diffusion coefficient [25]. 
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Another scenario, the Crushed Cloud model [26] invokes a "shocked" molecular cloud into 
which a radiative shock (typically v s ~ 100 km s _1 ) is driven by the high pressure of hot plasmas 
in the blastwave region. The cosmic -ray particles accelerated at a cloud shock are adiabatically 
compressed behind the shock front, resulting in enhanced synchrotron and 7T°-decay y-ray emis- 



sions Q27[] . It was shown by Ref. Q26H that reacceleration of pre-existing CR electrons and protons 
is capable of explaining both radio and GeV y-ray spectra using reasonable sets of parameters. The 
y-ray luminosity of Ly ~ 10 35 erg s _1 in 1-10 GeV agrees well with the theoretical expectation. 
In this scenario, the radio and y-ray emissions from the MC-SNRs are intimately connected, and 
therefore the possible radio-gamma correlation can be readily explained. 

The Runaway CR and Crushed Cloud models are not mutually exclusive. The former presumes 
unshocked molecular clouds as the y-ray production sites while the latter considers shocked clouds. 
Both unshocked and shocked clouds could be the important sites of the y-ray production. Also, in 
addition to the pre-existing GCRs in a molecular cloud, runaway CRs released from a shock may 
make a large contribution to the cosmic-ray density in the molecular cloud particularly at TeV 
energies. When such a cloud is struck and compressed by the blastwave, the y-ray flux at high 
energies should be enhanced with respect to what is presented in the model of Ref. [26]. 



Acknowledgments 

The Fermi LAT Collaboration acknowledges support from a number of agencies and institutes 
for both development and the operation of the LAT as well as scientific data analysis. These include 
NASA and DOE in the United States, CEA/Irfu and IN2P3/CNRS in France, ASI and INFN in Italy, 
MEXT, KEK, and JAXA in Japan, and the K. A. Wallenberg Foundation, the Swedish Research 
Council and the National Space Board in Sweden. Additional support from INAF in Italy and 
CNES in France for science analysis during the operations phase is also gratefully acknowledged. 



References 

[1] Blandford, R., & Eichler, D., Particle acceleration at astrophysical shocks: A theory of cosmic ray 
origin, Physics Reports, 154, 1-75, 1987. 

[2] Reynolds, S.P., Supernova Remnants at High Energy, ARAA, 46, 89-126, 2008. 

[3] Uchiyama, Y., Aharonian, F.A., Tanaka, T., et al., Extremely fast acceleration of cosmic rays in a 
supernova remnant, Nature, 449, 576-578, 2007. 

[4] Aharonian, FA., Akhperjanian, A.G., Aye, K.M., et al., High-energy particle acceleration in the shell 
of a supernova remnant, Nature, 432, 15-11 , 2004. 

[5] Esposito, J.A., Hunter, S.D., Kanbach, G., Sreekumar, P, EGRET Observations of Radio-Bright 
Supernova Remnants, ApJ, 461, 820-827, 1996. 

[6] Torres, D.F., Romero, G.E., Dame, T.M., et al., Supernova remnants and gamma-ray sources, Phys. 
Rep., 382, 303-380, 2003. 

[7] Abdo, A. A., Ackermann, M., Ajello, M., et al., The First Fermi Large Area Telescope Catalog of 
Gamma-ray Pulsars, ApJS, 187, 460-494, 2010. 

[8] Atwood, W.B., Abdo, A.A., Ackermann, M., et al., The Large Area Telescope on the Fermi 
Gamma-Ray Space Telescope Mission, ApJ, 697, 1071-1102,2009. 



6 



GeV Gamma-Rays from SNRs 



Yasunobu Uchiyama 



[9] Abdo, A. A., Ackermann, M., Ajello, M., et al., Fermi LAT Discovery of Extended Gamma-Ray 
Emission in the Direction of Supernova Remnant W51C, ApJL, 706, 1-6, 2009. 

[10] Abdo, A. A., Ackermann, M., Ajello, M., et al., Fermi-LAT Discovery of GeV Gamma-Ray Emission 
from the Young Supernova Remnant Cassiopeia A, ApJL, 710, 92-97, 2010. 

[11] Abdo, A. A., Ackermann, M., Ajello, M., et al., Gamma-Ray Emission from the Shell of Supernova 
Remnant W44 Revealed by the Fermi LAT, Science, 327, 1 103-1 106, 2010. 

[12] Abdo, A. A., Ackermann, M., Ajello, M., et al., Observation of Supernova Remnant IC 443 with the 
Fermi Large Area Telescope, ApJ, 712, 459-468, 2010. 

[13] Abdo, A. A., Ackermann, M., Ajello, M., et al., Fermi Large Area Telescope Observations of the 
Supernova Remnant W28 (G6.4-0.1), ApJ, 718, 348-356, 2010. 

[14] Abdo, A. A., Ackermann, M., Ajello, M., et al., Fermi-LAT Study of Gamma-ray Emission in the 
Direction of Supernova Remnant W49B, ApJ, 722, 1303-131 1, 2010. 

[15] Koo, B.-C, & Moon, D.-S. Interaction between the W51C Supernova Remnant and a Molecular 
Cloud. II. Discovery of Shocked CO and HCO+, ApJ, 485, 263-269, 1997. 

[16] Hewitt, J.W., Yusef-Zadeh, F., Wardle, M., et al., Green Bank Telescope Observations ofIC 443: The 
Nature of OH (1720 MHz) Masers and OH Absorption, ApJ, 652, 1288-1296, 2006. 

[17] Koo, B.-C, & Moon, D.-S., Interaction between the W51C Supernova Remnant and a Molecular 
Cloud. I. HI 21 Centimeter Line Observations, ApJ, 475, 194-210, 1997. 

[18] Green, D.A., A revised Galactic supernova remnant catalogue, Bulletin of the Astronomical Society 
of India, 37, 45-61,2009. 

[19] Castro, D., & Slane, P., Fermi Large Area Telescope Observations of Supernova Remnants Interacting 
with Molecular Clouds, ApJ, 717, 372-378, 2010. 

[20] Bykov, A.M., Chevalier, R.A., Ellison, D.C., & Uvarov, Y.A., Nonthermal Emission from a Supernova 
Remnant in a Molecular Cloud, ApJ, 538, 203-216, 2000. 

[21] Aharonian, F.A., & Atoyan, A.M., On the emissivity of n° -decay gamma radiation in the vicinity of 
accelerators of galactic cosmic rays, A&A, 309, 917-928, 1996. 

[22] Gabici, S., Aharonian, F.A., & Casanova, S., Broad-band non-thermal emission from molecular 

clouds illuminated by cosmic rays from nearby supernova remnants, MNRAS, 396, 1629-1639, 2009. 

[23] Ohira, Y., Murase, K., & Yamazaki, R., Gamma-rays from molecular clouds illuminated by cosmic 
rays escaping from interacting supernova remnants, MNRAS, 410, 1577-1582, 2011. 

[24] Aharonian, F, Akhperjanian, A.G, Bazer-Bachi, A.R., et al., Discovery of very high energy 

gamma-ray emission coincident with molecular clouds in the W 28 (G6. 4-0.1) field, A&A, 481, 
401-410, 2008. 

[25] Fujita, Y, Ohira, Y, & Takahara, F, Slow Diffusion of Cosmic Rays Around a Supernova Remnant, 
ApJL, 712, 153-156, 2010. 

[26] Uchiyama, Y, Blandford, R., Funk, S., Tajima, H., & Tanaka, T., Gamma-ray Emission from Crushed 
Clouds in Supernova Remnants, ApJL, 723, 122-126, 2010. 

[27] Blandford, R., & Cowie, L., Radio emission from supernova remnants in a cloudy interstellar 
medium, ApJ, 260, 625-634, 1982. 



7 



